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A set of stacked two-dimensional electron systems in a perpendicular magnetic field exhibits a 
three-dimensional version of the quantum Hall effect if interlayer tunneling is not too strong. When 
such a sample is in a quantum Hall plateau, the edge states of each layer combine to form a chiral 
metal at the sample surface. We study the interplay of interactions and disorder in transport 
properties of the chiral metal, in the regime of weak interlayer tunneling. Our starting point is a 
system without interlayer tunneling, in which the only excitations are harmonic collective modes: 
surface magnetoplasmons. Using bosonization and working perturbatively in the interlayer tunneling 
amplitude, we express transport properties in terms of the spectrum for these collective modes, 
treating electron-electron interactions and impurity scattering exactly. We calculte the conductivity 
as a function of temperature, finding that it increases with increasing temperature as observed 
in recent experiments. We also calculate the autocorrelation function of mesoscopic conductance 
fiuctuations induced by changes in a magnetic field component perpendicular to the sample surface, 
and its dependence on temperature. We show that conductance fluctuations are characterised by a 
dephasing length that varies inversely with temperature. 

PACS numbers: 73.20.-r, 73.23.-b, 72.20.-i, 73.21.Ac 



X 
S3 



I. INTRODUCTION 

Multilayer quantum Hall systems offer a setting in 
which to study the influence of electron-electron inter- 
actions and impurity scattering on tunneling between 
quantum Hall edge states. Specifically, consider a lay- 
ered conductor in a magnetic field that is perpendicular 
to the layers, with the field strength chosen so that a 
single layer in isolation would have quantised Hall con- 
ductance. Then, if interlayer tunneling is not too strong, 
the multilayer system exhibits a three-dimensional ver- 
sion of the quantum Hall effect and the bulk is insulat- 
ing at low temperatures. Under these conditions, edge 
states are present in each layer at the sample surface 
and are coupled by interlayer tunneling to form a surface 
phase, which is a chiral, two-dimensional metal,-'^'^ The 
contribution of this surface phase to the interlayer elec- 
tron transport properties of such systems has been iso- 
lated in experiments on semiconductor multilayers,^ and 
is dominant if samples are sufficiently small and cold. 

The consequences of impurity scattering for transport 
in the chiral metal have been discussed extensively from 
a theoretical viewpointitSdiSiSiLSiS and have been probed 
experimentally in several waysi^i-^^i-^-^i-^^i-^^i-^'^i-^^i-^^i^'^i-^^i-^^ 
Crucially, the chiral motion of electrons along the layer 
edges means that localisation is suppressed-iii^ As a re- 
sult, the surface conductivity in the interlayer direc- 
tion has a low-temperature limit that is non-zero, even 
though its measured value may be much smaller than 
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Separately, theoretical discussions of conduc- 



tance fluctationsiiSiSiLS have examined both their depen- 
dence on geometry in fully phase-coherent samples, and 
their dependence on the inelastic scattering length when 



this is smaller than sample size. Observations of re- 
producible mesoscopic conductance fluctuationsjiSii^ in- 
duced by small changes of magnetic field within a quan- 
tum Hall plateau, demonstrate that interlayer hopping is 
quantum-mechanically coherent and also provide a way 
to determine the inelastic scattering length. In addition, 
magnetoresistance in response to a field component per- 
pendicular to the sample surface has been proposed® and 
.|^gg(^ji4ji5ji7 ^ method for measuring the elastic scat- 
tering length. 

In contrast to these studies of disorder effects, past 
theoretical work on effects due to electron-electron inter- 
actions in the chiral metal has been limited. There have 
been discussions, first, of the temperature dependence 
of the inelastic scattering length-'- and, second, of the 
fact that there is no zero-bias anomaly in the tunneling 
density of states (or any related contribution to the con- 
ductivity), because of ballistic motion of charge in the 
in-layer directioniSiS 

Against this background, recent experiments find- 
ing a significant temperature dependence to the surface 
conductivitjiiSii^ are striking as likely indications of in- 
teraction effects, and provide one of the motivations for 
the work we present here. In particular, the fact that 
conductivity is observed to increase with increasing tem- 
perature presents a puzzle for theory. Some straightfor- 
ward potential explanations are specifically excluded by 
the experimental design: large ratios of sample perime- 
ter to cross-sectional area ensure that surface states make 
the dominant contribution to the measured conductance; 
and sample perimeters much longer than the inelastic 
scattering length ensure that weak localisation effects are 
absent. For samples studied in Ref. the measured 
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conductivity cr(T) increases by about 7% in the temper- 
ature range from 50mK to 300mK, implying a tempera- 
ture scale of cr(T) • [da{T) / dT]-^ - 4K, which is similar 
to that for other interaction effects in quantum Hall sys- 
tems 

In this paper we study interactions and disorder in the 
chiral metal, working in the experimentally-relevant limit 
of weak interlayer tunneling. Treating tunneling pertur- 
batively, Coulomb interactions and impurity scattering 
can be handled exactly by means of a straighforward ap- 
plication of bosonization. We calculate the full tempera- 
ture dependence of the conductivity. We also study con- 
ductance fluctuations induced by magnetic field changes, 
obtaining their autocorrelation function and its depen- 
dence on temperature. Making appropriate parameter 
choices, our results for both quantities are consistent with 
experimental findings. A short account of this work has 
been presented previously, in Ref. 12(1 

Our work differs from most of the extensive literature 
on tunneling between quantum Hall edges states in two 
important ways. First, while much previous work has 
been concerned with edge states of fractional quantum 
Hall systems (SiiSSiSiSi including multilayer s&wples^^^ 
our focus is on the integer quantum Hall effect. Sec- 
ond, whereas most past work (with some exceptions: see 
Refs. 27, 28 29 30j has been restricted to systems with 
only short-range interactions, we find that the long-range 
nature of Coulomb interactions, which we treat in full, is 
central for the results we obtain. 

The remainder of this paper is organised as follows. We 
develop a model for the chiral metal in Sec. ITTIand show 
how bosonization can be used to give an exact description 
of the collective excitations. Sec. lHII contains calculations 
of the temperature dependence of the conductivity. We 
study conductance fluctuations in Sec. IIVI and discuss 
our results in Sec. 



II. MODELLING THE CHIRAL METAL 

In this section we summarise the physical ingredients 
that are important for modelling transport between edge 
states in multilayer conductors and set out the length- 
scales that characterise the system. We introduce a 
Hamiltionian in terms of fermionic operators for edge 
electrons. We bosonize this Hamiltonian, obtaining a 
result which is quadratic in boson operators if interlayer 
tunneling is omitted. Finally, we express the two-electron 
correlation function that is central to transport calcula- 
tions in terms of boson correlators. 



A. Ingredients, lengthscales, and parameters 

A multilayer conductor is illustrated in Fig^ We use 
coordinates with the cc-axis parallel to the layer edges, 
and treat a sample of N layers with layer index n and 
layer spacing a. Consider the system in the presence of 



a perpendicular magnetic field of strength B, with the 
chemical potential lying between the lowest and first ex- 
cited Landau levels. In the bulk of the sample single 
particle states at energies close to the chemical potential 
are localised by disorder. At the sample surface in this 
energy range, edge states propagate in the confining po- 
tential Vcdgc{y) at a velocity v. Interactions modify the 
confining potential and the edge velocity: we denote by 
vf the velocity allowing for Hartree contributions. Edge 
states have a width w in the y-direction, which is set by 
the magnetic length Zb in a clean sample, and by the 
bulk localisation length ^ in the presence of impurities. 
We use a one-dimensional decription of the edge state in 
each layer, projected onto the x-coordinate in the stan- 
dard way. 

Out theoretical treatment takes account only of one 
edge state in each layer and is therefore appropriate for 
a system in which electrons are spin polarised. In fact, 
some of the experiments we refer to, including those on 
the temperature-dependence of conductivity^ are for 
systems with Landau level filling factor per layer of = 2. 
It is appropriate to apply our theory to these systems pro- 
vided electrons with opposite spin directions contribute 
additively and incoherently to the conductivity. 

The system of edge states can be characterised using 
three lengthscales. First, impurities, which generate only 
forward scattering with a phase shift, result in an elas- 
tic mean free path l^h the distance over which a phase 
shift of order 27r is accumulated. Second, temperature 
T in combination with the velocity v-p can be expressed 
in terms of the thermal length Lt = hvp/k-QT. Third, 
interlayer tunneling with amplitude t± can be parame- 
terised by the characteristic distance lj_ through which 
electrons move in the chiral direction between tunneling 
events. The value of l± can be expressed in terms of the 
interlayer diffusion constant D: since, for small t±, in- 
terlayer hops are of length a and occur at a rate vf/1±, 
one has /_l = a^vp/D. In turn, this can be expressed 
in terms of the conductivity, using the Einstein relation 
and the fact that the density of states is n = X/inahv-p^ 
giving l± = a{e'^ /2nha).^ 

Parameter values for the experiments of Refs. 0, UtI 
and are as follows. Samples consist of A'' ~ 50 - 100 
layers with spacing a — 30nm. The mean free path is 
estimated^^ to be ^oi ^ 30nm. An upper bound on v-p, 
reached in samples with a steep confining potential is 
vp ~ wc^B, where uq is the cyclotron frequency. It has 
the value ujqIb = 1-7 x lO'^ms"! in GaAs at 6.75 T. With 
this value, ix ^ 10/im at T = lOOmK. Finally, for a 
surface conductivity of cr = 1.3 x lG^'^e^/27rfi, (which lies 
within the observed range aX v — 2), Ij^ = 40/im. We 
are therefore concerned with the regime lc\ <C L'p <^ l±, 
and this motivates our approach, based on a perturbative 
treatment of tunneling. 
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FIG. 1: A multilayer conductor, showing the orientation of 
axes in our coordinate system, with edge states propagat- 
ing in the a;-direction. The form of the confining potential 
Vedgeiy) is illustrated top left. Interlayer tunneling amplitude 
and spacing are denoted by t± and a, respectively. 



B. Fermionic Hamiltonian 

Our model Hamiltonian, Ti. — Tio + ^dis + ^hop + ^int , 
has single-particle terms Ho, Hdis and Tihop, represent- 
ing, respectively, free motion along each edge, impurity 
scattering and interlayer hopping, and a contribution 
Tiint from Coulomb interactions. We write it in terms 
of the electron creation operator cj„ for an edge state 
with wavevector q in layer n, taking sample perimeter L 
so that q = 27rng/L, where Uq is integer. The creation 
operator at a point is 



-iqx J 
^ '-qn ■ 



(2.1) 



We take Vn{x) to be Gaussian distributed with zero- 
range correlations and strength A: [Ki(a;)]av = and 
\yn{x)Vn' {x')]a.-v = ASn,n'^{x — x'). This disordcr term 
can be removed by means of a gauge transformation on 
the fermionic field operators, under which 



where 



Ux) 



hv 



dx'Vn{x') 



(2.6) 



(2.7) 



is the phase shift acquired under forward scattering from 
the impurities. The clastic scattering length is related 
to the disorder strength A by lc\ — H'^v'^/A. Under this 
gauge transformation. Ho + Hdis Hq. The hopping 
term, however, picks up a dependence on the disorder, 
and after the transformation is 

^hop dx[t^{n,x)^pi+^{x)Mx) +R. c], (2.8) 



where 



(2.9) 



We ignore the effects of this gauge transformation on the 
boundary conditions applying to t/^nix), which is justi- 
fied at temperatures large compared to the single-particle 
level spacing. With this, Hq + Hint is unaffected by the 
gauge transformation, and gauge transformed operators 
cj„ can be defined by inverting Eq. H2.1|l . All further 
references in this paper to fermionic operators are to the 
gauge-transformed ones. 



We normal order the Hamiltonian with respect to a vac- 
uum in which states are occupied for q < and empty 
otherwise. Then 



Ho — 



ihv^ J dx ■.^pl{x)dxtlJn{x): , (2.2) 



and 



Hhop = ^ j dx[ti_ijj\^^{x)i;n{x) + H. c] 



(2.3) 



The interaction contribution, written in terms of the pro- 
jected density p{x) — ipl^{x)ipn{x) with a two-particle po- 
tential Un-m(x — x'), is 

^int = j '^^ j dx' : Pn{x)Un^rn{x - x')pm{x') : . 

nm 

(2.4) 

Finally, writing the impurity potential projected onto 
the edge coordinate in the nth layer as Vn{x), we have 



^dis = X! / dxVn{x) : ip'l{x)ij;n{x) ■■ 



(2.5) 



C. Bosonised Hamiltonian 

We bosonize the Hamiltonian in the standard way, ex- 
pressing Ho + Hint in terms of non-interacting collective 
modes. Since Tihop transforms into a cosine function of 
the boson creation and annihilation operators, we treat 
it perturbatively. To justify this, we require that t_L 
should be small. Since t± is a relevant perturbation,^^ 
we also require that temperature should not be too small: 

Boson creation operators are defined in the usual way 
(see, for example, Ref. l3ll) as 



fot = I V c'^ c 

qin 1^ \i/2 / . '^r+q,m'^r,m 



(2.10) 



for g > 0. Fourier transforming the interaction potential 
and expressing the result as a velocity, we introduce 



Un-m{q) = {2Trhy 



dxe""'Ur, 



(2.11) 



The Fermi velocity renormalised by Hartree interactions 
is vf = V — u„(0), where the divergence which arises 
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in the sum in the case of Coulomb interactions is can- 
celled by contributions to v from a neutralising back- 
ground. The Hamiltonian in the absence of hopping (and 
omitting fermion number terms which appear at electron 
densities different from that of our vacuum) is 

T-io + Hint = X! X! ^[^F + Un~rn{q)]qblnbq„i ■ (2.12) 
mn g>0 

The combination Ho + Hint is diagonalised by Fourier 
transform in the layer index n. We impose periodic 
boundary conditions on n, define the wavevector k = 
2nkTT/Na, with integer and —n/a < k < n/a, and set 



For isolated layers, taking the limit of large a, the sum 
on p may be replaced with an integral and one recov- 
ers the dispersion relation of edge magnetoplasmons in a 
single layer system, known from previous work^SiS 

For the multilayer system the expression for the dis- 
persion relation may be simplified in two stages. First, if 
the layer spacing is small (a ^ w) the sum on p may be 
omitted, so that 



uj{q,k) = v-pq 1 + 



(2.20) 



If, in addition, interactions are weak (uj <C k ^) 



and 



qk 



1 ^ 
^ n— 1 



4g,fc) = ^e™'='^7.„(g). 



(2.13) 



(2.14) 



}{q,k) = upg 1 



(2.21) 



In the following we obtain detailed results for sys- 
tems with wide edges using the dispersion relation of 
Eq. (|2.20() . and for systems with narrow edges using the 
dispersion relation of Eq. H2.21|l . 



Then 

Ho+7^i„t-^^fia;(g,fc)6j,6^, (2.15) 

k q>0 

where the excitation frequencies are 

iuiq,k) = [vF+u{q,k)]q. (2.16) 

The Coulomb interaction, regularised at short dis- 
tances by a finite width w for edge states, has the form 



Unix) 



(2.17) 



The edge state width w is set by the localisation length ^ 
of localised states in the bulk of the sample at the Fermi 
energy. In a clean sample with well-separated Landau 
levels, ^ ^ ^B, but in a highly disordered sample with 
Landau levels that are broad in energy one may have 
^ ^ Zb- The value of w proves important in matching 
our results to experiment, as we discuss in Sec. 11110 

We write the Fourier transform, using the Poisson sum- 
mation formula, as 



u{q,k) =VF-^Y\ 



dxdz- 



^~i{qx+kz+2'iTpz /a) 
y/x^ + + 



(2.18) 



and find 



Loiq, fe) = vpg I 1 + K ^ Qp^e-"""^- 



(2.19) 



with Qp — q'^ + (k + 27rp/a)^ and p integer. Here, the 
inverse screening length k = /Airereohvpa characterises 
the interaction strength. 



D. Two-particle correlation function 

A central quantity in our calculations of transport 
properties is the two-fermion correlation function 

G{x,t) ^ (V't(a;,i)V'„+i(x,i)^l+i(0,0)V^„(0,0)), 

(2.22) 

where (...) = Tr(e^'^^ . . .)/Tr(e~^'^) and operators are 
written in the Heisenberg representation, with 0{t) = 
^iHt/hQ^-iHt/h^ evaluate this in the absence of tun- 
neling, so that H = Ho + Hint- 

As a first step, define the boson field operator— 

Mx) = -J2 i^-'^'^u + e'''"6,„) (2.23) 

where e is a short-distance cut-off. Omitting Klein factors 
(which cancel from G{x,t)), the fermion and boson field 
operators are related by 

M^) = (27re)-i/2 {-i^nix)) . (2.24) 
The correlation function is 
G{x,t)-- ^ 



(27re)2 

We define its logarithm S via 
G{x,t) 



gi0„(x,t)g-i0„ + i(x,t)gi0„ + i(O,O)g-i0„(O,O)\ 



(2.25) 
(2.26) 



Because H is harmonic, S can be expressed as 

S = -^{{Mx,t)-^n+l{x,t) + (f>n+l{0,0)-(pn{0,0)f) 

+ i [0„ {x,t)- </.„+! {x,t), 0„ (0, 0) - 0„+i (0, 0)] 
- 2 log e. 

(2.27) 
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The thermal average and the commutator appearing 
in this expression can be evaluated in the standard 
way via a mode expansion, by expressing 0„(x,t) in 
terms of boson creation and annihilation operators using 
Eq. (|2.23|1 . Taking the thermodynamic limit and so re- 
placing wavevector sums with integrals, with f3 = \/k^T, 
we arrive at 

S{x,t,T) = -2\og€- - dk{l- cos ak) -^e'"^ 

T^J-Tr/o Jo Q 

X ( coth{ l3 huj {q, k) / 2 cos {qx~ uj{q,k)t)] (2.28) 



+ i sin {qx— uj{q, k)t) 

It is useful to note that 

G{~x,-t)^G{x,t)*, (2.29) 
and also to define a frequency-dependent correlator, 

G{x,n) = [ dte'"'G{x,t). (2.30) 



III. CONDUCTIVITY 

In this section we express the conductivity cr(T) ob- 
tained from a Kubo formula in terms of the two-fermion 
correlation function calculated in Sec. IIIDI We also set 
out the steps required for a numerical evaluation of cr(T), 
present our results, and compare them with the experi- 
mental data of Ref. llSL 



A. Kubo formula for conductivity 

The operator for the interlayer current density between 
layers n and n -I- 1 is 

jn{x) = ^ (^tj_{n, x)^l_^_^{x)ipn{x) - H. c.) . (3.1) 

The real part of the conductivity at frequency Q is given 
by the Kubo formula — 

a{n, T) = ^1 sin nt f dx f dx' 

X (7„(x,t)j„(x',0)) . (3.2) 

To leading order, the interlayer hopping appears only in 
the current operators, and we evaluate the thermal aver- 
age using a Hamiltonian from which interlayer hopping 
is omitted. 

Substituting for j„(a;, t) using Eq. H3.1|) gives an expres- 
sion for the conductivity of the chiral metal with a given 



configuration of disorder: to leading order in tj_(n, x), 
X t±(n, x)t*j^(n, x) 

X {Tpl^{x, 1)11) n+l{x, i)V'l+i(x',0)^„(a:',0)). 

(3.3) 

Averaging over disorder configurations yields 

[ij_(n,a:)il(n,x')]av = tie-l"l/'"' (3.4) 

and hence 



a(r!,T) 



dx 



dt sin U,t 



h nh? J 2?ci 
X (V'i(x,t)V„+i(x,t)Vl+i(0,0)V„(0,0)). (3.5) 



This result can be expressed in terms of the time or fre- 
quency dependent two-particle correlation functions de- 
fined in Sec. Ill D1 Setting = we find 



8nalc\t'^f f dx 



ft2 



2L 



Re 



dnG{x,n)\^^^ 



For a boson dispersion relation ^(q, k) — vpq, as results 
from the Hartree approximation, the fermion correlation 
function factorises into independent contributions from 
each layer. These have the form 



(^t(x,t)V„(0,0)) = ^ 



dk- 



l + e 



f3hvpk 



(3.7) 



and we find a temperature-independent conductivity 



ain,T) 



h h^vi i + my 



2;2 /„2 ' 



(3.8) 



cl/ ""F 



which in the zero-frequency limit has the value 

(3.9) 



2t\ leia 



h h^vl 



More generally, with an arbitrary boson dispersion re- 
lation a simplification of Eq. (|3.t)|l is possible for 4i «C Ltj 
since G{x, t) varies with x only on the scale Lt while the 
correlator [tj^(n, a:)t^(n, a;')]av has range ^oi- We get 



/C30 
dttl\aG{Q,t) 
-OO 



(3.10) 



B. Evaluation of a{T) 

To find the temperature dependence of the conductiv- 
ity we must combine Eqs. (|?!^ . and ifXTUI) . A 
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first step before numerical evaluation is to isolate the de- 
pendence on the cut-off e and take the limit e — > 0, as we 
describe in this subsection. 

We start from the expression given in Eq. (|2.28|) for the 
logarithm of the two-particle correlation function, which 
we evaluate at a; = 0. It is convenient to separate out a 
zero-temperature contribution by writing 



S{t,T) = S{t,0) + AS{t,T) 



(3.11) 



and also to split S{t, 0) into real and imaginary parts, 
with 

Sit,0) = U{t)-iV{t), (3.12) 
where h{{t) and V{t) are real for t real. Then, writing 



a{T)=a{0) + AaiT), 
we obtain from Eq. (|3.1U|) 



a{0) 







(3.13) 



(3.14) 



and 



Aa(T)= 



dUe"(*)sinV(t)[e^^(*'^)-ll . (3.15) 



In the case of a linear boson dispersion relation, uj{q, k) 
vfQ, the functions U{t) and V{t) have the forms 



Viin(t) = vr- 2tan"^(e/i;Ft). 



(3.16) 
(3.17) 



Adding and subtracting these expressions from the ones 
for U{t) and V{t) with a general dispersion relation, we 
find 

a r'"' 

U{t) ^Uiin{t) + - dk{l-cosak) 

J-TT/a 



/ — e^'^'^[cos {lo {q, k)t) ~ cos (vFqt)] 
Jo 1 



(3.18) 



and 



V{t) = Vii„(<) + - / dk{l - cosafc) 

J-TT/a 

X / — e^*^^ [sin {Lu{q, k)t) — sin {v^qt)] 
Jo 1 

Finally, we have 

AS{t,T) = / dk(l- cos ak) 

TT J-7T/a 



(3.19) 



(3.20) 



dq 



X / — e "'^{1 - cos to {q,k)t) 
la 1 



coth 



^ Pf^{q,k) 



The advantage of casting the equations for the con- 
ductivity in this form is that the momentum integrals in 



Eqs. (|XTH|l . l|TT^ and can be performed at e = 0, 

since the integrands decay fast enough at large q for con- 
vergence. Dependence on e is confined for small e to the 
functions U\in{t) and V\in{t), and from Eqs. H3.16|l and 
(|3.17l) one sees that it is important only for t ~ 0(e). It 
is therefore convenient to separate the integration range 
in Eq. (|3.14|) into two parts, <t < R and R<t<oo, 
with e ^ i? <C 1. In the first interval U{t) = U\in{t) and 
V = Viin(i); in the second interval one can set e = 0. 

Let the contributions to a(0) from the two intervals be 
cr^^^ and a'^^\ Writing t' — vpt/e we have 



(1) _ 2e Q-Q 



dt't' e^'-''hmV{t') (3.21) 



which gives 



.(1) 



2fTo 



2t' 



dt't' - 

TT ./o 1 + 1 + 



CTq . 



(3.22) 



Evaluation of cr'^-* requires a numerical calculation, and 
we present results in Sec. IIII Dl 

Finally, turning to the conductivity at non-zero tem- 
perature, we note that there are no extra difficulties in 
the evaluation of Act using Eq. (|3.15|l . The function 
AS'(f,T), can be computed numerically with e = 0, and 
AS{t,T) ^ as i 0, so that ACT(r) has no contribu- 
tion from the integration interval < t < R in the limit 
e ^ 0. 

In summary, when evaluating cr(0) or Aa{T) using 
Eqs. ifXTljl and l|XTH|l . the functions l({t), V{t), and 
AS{t, T) may be evaluated numerically by setting e = 
in Eqs. H3.18|l . I|3.19|l . and H3.2()(l . and the results used in 
Eq. (jXT^ to find a^'^\ To this one must add a^'^^ = in 
order to obtain the zero temperature conductivity a{Q). 
These equations combine with Eq. (|3.15() for Aa{T) to 
give a computationally tractable, though non-trivial, ex- 
pression for cr(r). 



C. Conductivity at zero temperature 

The conductivity at zero temperature and zero fre- 
quency is determined solely by the low energy limit of the 
group velocity for excitations, since no other modes are 
excited as T, 17 ^ 0. This zero frequency limit is reached 
as q, the wavevector component in the chiral direction, 
approaches zero. The group velocity, duj{q,k)/dq\q^o = 
VFa{k), is in general a function of fc, the wavevector com- 
ponent in the interlayer direction. 

To determine (t(0), a useful procedure is to consider 
a model dispersion relation which is exactly linear in 
q: io{q,k) = VFqa(k). A linear dispersion relation is 
also of interest in its own right. It arises from an in- 
teraction that in real space is short range in the chi- 
ral direction, x: Un{x) = gnS{x), giving a{k) = 1-1- 
{2TThvF)~^ e*'^"'^5„. With a linear dispersion relation, 
q-integrals in the expressions leading to G(x, t) can be 
evaluated analytically, greatly simplifying the calculation 
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of conductivity. As we show in the following, for the limit 
lei -^T that we consider, a dispersion relation linear in 
q yields a temperature- independent value of conductivity. 
For interactions, such as the Coulomb potential, that are 
not short range in x, linearisation of the dispersion rela- 
tion gives only an approximation to G{x, t). The value of 
(t(0) that results from integrating this approximate form 
for G{x, t) is nevertheless exact (at the leading order in 
t± considered throughout this paper). This fact is clear 
on physical grounds, since we have correctly accounted 
for the dispersion relation at low energy. It may also be 
derived formally, as follows. 

Starting from Eq. (jXTUI) . we deform the contour for the 
time integral into the semicircle at infinity in the lower 
half of the complex plane, writing t — Ir + iti with tR 
and tj real. Then in Eq. (|2.28|l we have the factor 

dq- exp(— — iqx — itRUj{q, k) -f tjui{q, k)) . (3.23) 
1 

This must be evaluated for all values of t lying on the 
deformed time integration contour. When \tfi\ is large, 
exp(— i<jjw(q, k)) is a rapidly oscillating function of q, and 
the g-integral can be computed using the method of sta- 
tionary phase: since io{q, k) is a monotonically increasing 
function of g, the dominant contribution comes from the 
vicinity of the end-point at g = 0. Similarly, when </ is 
large and negative, eyi'p{tiuj{q, k)) is small for most values 
of q, and the g-integral can be computed using steepest 
descents: again, the dominant contribution comes from 
the vicinity of g = 0. In both instances we may approx- 
imate uj{q,k) by its form linearised about q — 0; after 
linearisation the g-integral can be evaluated analytically. 
This calculation yields 



G(0,i) 



1 



1 



1 



(27r)2 Vsinh {TTt/(]h) J w|(e -|- itf 



(3.24) 



X exp / dk{\ — cos ak) log 

\ ^ J, 

Substituting this into Eq. H3.10|l we obtain 

2(T / 2{i z*^/^ \ 

a(T) — exp / dfc (1 — cosa/c) loga(fc) 

\ Ja J 

f dtee / nt/(3H V 
J (e2 + i2)2 {sinhint/ph)) ' 

In the limit e ^ 0, the t integral gives tt/2 regardless 
of temperature, demonstrating that, for systems with a 
linear dispersion relation, in the regime Id ^ ^Ti cr(T) 
is independent of T. We find 

(2fl z*^/^ \ 
/ (ifc(l — cosafc) loga(A:) 1 . 
^ Jo J 

(3.26) 

This is our final result for the dependence of a{Q) on the 
dispersion relation as parameterised by a{k). 



D. Results 

We are now in a position to calculate the conductiv- 
ity for a system with Coulomb interactions by evaluating 
numerically the formulae we have derived: first, the zero- 
temperature value using the results from Sec. IIII 01 and 
then the full temperature-dependent conductivity using 
the results from Sec. IIII Bl We investigate variation of 
the conductivity with two parameters, the Fermi veloc- 
ity v-p and the edge state depth w, and seek values of 
these parameters for which our results match the exper- 
imental data of Ref. 18. The parameters enter the dis- 
persion relation w(g, fc) directly, and v-p also appears in 
the inverse screening length k. The interaction strength 
is set by the combination Ka (recall that a is the layer 
spacing). A scale for temperature is set by vp and a, 
via To = hvp/akB, so that T/Tq ~ q/Lt- A scale for 
conductivity is given by ao, its value in the Hartree ap- 
proximation. 

At a qualitative level, the effect of interactions on the 
conductivity can be anticipated by starting from the ex- 
pression given in Eq. H3.9|l for this quantity within the 
Hartree approximation. In turn, that expression can be 
understood in terms of a calculation of the interlayer tun- 
neling rate, based on the Fermi golden rule: the rate 
involves the square of a matrix element between initial 
and final states on adjacent layers, and a power of the 
density of states for both the initial and the final states. 
The squared matrix element, allowing for disorder which 
affects phases of initial and final states separately, con- 
tributes a factor of t'^ld to ctq- The form of the density of 
states on a single edge, 1/2'Khvp, implies that (Tq oc Wp^. 
Returning to a full treatment of the interacting system, 
we note that the effect of interactions is to generate an 
energy-dependent group velocity in place of a constant 
value, vp. In effect, the value of cr(T) at a particular 
temperature involves a thermal average of the inverse 
square of the group velocity. Because Coulomb inter- 
actions increase the group velocity at low energy, they 
decrease conductivity at low temperature; equally, be- 
cause the group velocity approaches vp at high energy, 
the conductivity approaches ctq at high temperature. 

Turning to detailed results, the dependence of cr(0) 
on w/a and Ka is shown in Fig. |21 as obtained from 
Eq. H3.26|) using a{k) = 1 -|- Ke"™''^'/!^!- Interactions 
reduce the value of the conductivity, by a factor which 
is large if Ka is large. The variation of cr(T) with T is 
illustrated in Fig. El for a system with the dispersion re- 
lation appropriate for narrow edge states, Eq. (|2.21|l . In 
this case the k integrals in Eqs. (|3.18l) and H3.19|l can 
be done analytically, leaving only the g and t integrals 
to be evaluated numerically. Finally, the behaviour of 
(t(T) for a system with wide edge states {w > a) is pre- 
sented in Fig. ^ In this case the dispersion relation is 
as given in Eq. H2.20|l . analytical progress does not seem 
possible, and integrals on fc, g and t must be evaluated 
numerically to obtain cr(T). We note in passing that we 
checked that there are only small changes to the results 
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FIG. 2: Conductivity at zero temperature, as a function 
of interaction strength, parameterised by inverse screening 
length K., for various edge state widths w. 



1 




10 20 30 



FIG. 3: Dependence of conductivity on temperature for nar- 
row edge states, with interaction strengths na = 1 and na — 5. 




FIG. 4: Dependence of conductivity on temperature for wide 
edge states with w = 4a and interaction strengths na — 1 and 
na — 50. 



presented when using the more complete form of the in- 
teraction given in Eq. I|2.19|l . including the sum on p. 

Examining these results, it is evident that the gen- 
eral shape of (t{T) does not vary greatly with param- 
eters: the temperature dependence is quadratic at low 
temperatures, has a roughly linear region at intermedi- 
ate temperatures, and approaches ctq in the high temper- 
ature limit. The quadratic dependence at low temper- 
ature is universal, but the extent of the roughly linear 
region at intermediate temperature is model-dependent. 
Moreover, scales in this temperature dependence change 
dramatically with parameter values. The value of the di- 
mensionless temperature T/Tq at the crossover between 
the low and intermediate temperature regimes is depen- 
dent on K (see Fig. ^ and varies even more strongly with 
w (compare Figs. |21and^. In addition, the magnitude 
of the variation in cr(T) between low and high T depends 
very much on the values of w and Ka. In order to repro- 
duce the experimental observation of a nearly linear in- 
crease in (t(T), by about 7% between the temperatures of 
50mK and 300mK,-lS' we require parameters which place 
the experimental temperature window in the intermedi- 
ate regime for behaviour, so that quadratic variation of 
cr(T) with T occurs only in a temperature range below 50 
mK, and saturation of cr(T') occurs only above 300mK. 
Since the available data is not sufficiently detailed to jus- 
tify a formal fitting procedure, we instead survey the con- 
sequences of a range of parameter choices in our results 
and examine the match to experimental observations. 

We begin by considering narrow edges states, using 
the results shown in Fig. |3| Supposing vf ~ ujqIb, 
which represents an upper bound on vf, we have vf — 
1.7 X lO^ms"^. With a ~ 30nm, we find ko ~ 1 and 
To ~ 40K. Taking these values, the variation in cr(T) 
over the experimental temperature range is very small 
and quadratic, in disagreement with observations. A re- 
duction in the value of wp serves to decrease the tem- 
perature scale To, and also increases k. It is possible to 
generate approximately linear variation of a{T) with T 
in the experimental temperature range by using a suffi- 
ciently small value of vf (reduced from the upper bound 
by ^ 0(10^)), but we know of no reason for vp to be so 
small. 

We therefore turn to theoretical results for wide edge 
states, as illustrated in Fig. 0] In this case, we find that 
large values of w greatly reduce the temperature range 
over which a{T) varies quadratically with T, and can 
lead to approximately linear variation in the experimen- 
tal temperature range. A second consequence of large 
w is that the conductivity change cr(oo) — cr(0) is re- 
duced. This tendency can be counteracted by increasing 
the interaction strength Ka. We find that observed be- 
haviour can be reproduced by taking w = Aa — 120nm 
and Dp = 3 X lO'^ms""'^ (giving na = 50). The tempera- 
ture dependence of cf(T) obtained using these parameter 
values is shown in Fig.jSjfor temperatures below 400mK. 

This choice of parameters, and its implications, merit 
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FIG. 5: Dependence of conductivity on temperature for 
w — 4a, with a = 30nm, i;f = 3 x lO^ms"^ and ctq ~ 
1.893 X 10~^e^/27rft (full line), compared with experimental 
data (points) taken from Fig 2 of Ref. (data set for Fractal 
2). 



further discussion. First, we note that there are two sep- 
arate experimental indications that edge states have a 
width closer to the value we have adopted, of 120nm, 
than to the conventionally expected value of — lOnm. 
One comes from measurements of bulk hopping transport 
in multilayer samples"^^, which give a localisation length 
of ^ = 120nm: one expects w ~ f . The other comes 
from studies of conductance fluctuations discussed in 
Sec. IIVI These yield a value for the inelastic scattering 
length, from the amplitude of fluctuations, and a value 
for the area of a phase-coherent region perpendicular to 
the applied field, from the correlation field for fluctua- 
tions. The ratio of this phase-coherent area to the inelas- 
tic scattering length implies an edge state width which 
is also much larger than w ~ 70nm. Next, turning 
to the value of v-p , which we have taken 50 times smaller 
than for edge states in a steep confining potential, we 
note that large edge state width favours a small value for 
vp, because wide edge states penetrate into the bulk of 
the sample where both the confining potential gradient 
and the drift velocity of electrons moving in this poten- 
tial are small. Finally, we comment on the fact that 
accepting a small value for v-p implies a large value for 
ctq, if other parameters are unchanged. In fact, large w 
acts in the opposite direction, to reduce the effective tun- 
neling amplitude t± between edge states, since different 
portions of the edge contribute to the amplitude with 
different phases, so that there are partial cancellations. 
To account for the magnitude of the measured^^ conduc- 
tivity, 1.5 X 10~^e^/27rfi,, using the value for the mean 
free path 1^ = 30nm derived from magnetoresistance 
measurements^^ requires an effective value of t± about 
50 times smaller than bare estimate^, of 0.12 meV. This 
is a surprisingly strong supression of tunneling, though 
possible if edge states in successive layers have different 



displacements from the surface, as suggested in Ref. 

IV. CONDUCTANCE FLUCTUATIONS 

It is found experimentally that mesoscopic fluctua- 
tions in the conductance of the chiral metal are induced 
by small changes of magnetic field within a quantum 
Hall plateauiiSiiS These conductance fluctuations are ob- 
served in samples with a perimeter that is several times 
larger than the estimated inelastic scattering length. Un- 
der such conditions, it is not initially clear why the mag- 
netic field component perpendicular to layers in the sam- 
ple should influence conductance in this way, since in the 
simplest picture electron trajectories enclose flux only by 
encircling the sample. More realistically, a number of 
possibilities are evident the sample walls may lie at an 
angle to the layer normal, either on average or because 
of surface roughness, or finite edge state width may be 
important. In our theoretical treatment of conductance 
fluctuations we avoid specific assumptions about this as- 
pect of the system by considering fluctuations that result 
from variations in a magnetic field component B± per- 
pendicular to the sample surface. The amplitude of fluc- 
tuations is not affected by this choice. By contrast, the 
scale for the correlation field of fluctuations is dependent 
on the model chosen for flux linkage. 

In a general setting, there are two possible reasons for 
the amplitude of conductance fluctuations to decrease 
with inceasing temperature. One is because of a decrease 
in the inelastic scattering length; the other is because of 
thermal smearing. In the case of a chiral metal only the 
flrst mechanism operates, because states at different ener- 
gies are perfectly correlatedi— In this sense, conductance 
fluctuations offer a rather direct probe of interaction ef- 
fects. 

In this section, in place of conductivity a, we are con- 
cerned with the conductance g = aL /Na of a finite sam- 
ple and fluctuations 6g — g— [g]av about its average value. 
We denote the average within the Hartree approximation 
by go = aoL/Na. We derive an analytic expression for 
the autocorrelation function of conductance fluctuations 
induced by B±. We focus on its temperature depen- 
dence at low temperatures, obtaining a scaling form for 
the regime in which cr(T) « cr(0). We compute the scal- 
ing function, evaluate our expressions numerical ly, and 
compare our results with the observations of Ref. Il9i 

A. Correlation function 

The conductance autocorrelation function 

Fi6B) = [5g{B^)5g{B^ + 5B)]^, (4.1) 

is characterised by the amplitude F{Q) and by the corre- 
lation field. An obvious field scale is set by a flux density 
of one flux quantum <f>o through a rectangle with sides 
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proportional to the layer spacing and the thermal length, 
and we define Bq = ^o/SttoLt — fi/eaLT^. We also in- 
troduce a dimensionless field variation b = SB/Bq, which 
depends on temperature through L^, and a temperature- 
independent reduced field h which has dimensions of 
wavevector: h — 6/Lt = e5B/ah. 

With a suitable choice of gauge, the transverse field 
enters the Hamiltonian only as a phase for interlayer hop- 
ping. Taking for convenience B± = 0, in the presence of 
non-zero SB Eq. 1)2. 9|l is modified to 



(4.2) 



This additional, field-dependent phase alters Tihop and 
consequently the current operator. 

An expression for the conductance of a sample with 
a specific disorder configuration is obtained by scaling 
Eq. (|3.3|) with the sample dimensions. Taking account of 
the field-dependent phases in the current operator and 
substituting into the definition of F{SB), after some ma- 
nipulation we arrive at 

nSB) = fd^' f dy' (4.3) 



X J dtitG{x-x',t) J dt' it' G{y ~ y\t') 

X ^g«'»(^-^') _|_ g-fi(2;-2;')~j ^C(x,x') ^C{y,y') 

Two contributions to this expression arise from the dis- 
order average: 



C{X,X')^- i [(^?„+l [X) - On (X) - On+l {x')+en{x)f] 



(4.4) 



and 



D„m {x, x'; y, y') = [{9n+i{x) - On [x) - 9n+i{x') + 9n [x')) 
x((?„,+i(y)-0,„(y)-0,„+i(y') + e™(2/'))]av- (4-5) 
Both may be evaluated using the result (for x, y > 0) 



iix)9miy)]a,v = -j^ mm{x, y} 



The equation for C gives 



^C(x,x') _ ^~\x~x'\/l^i 



(4.6) 



(4.7) 



which in the limit of small ^ci can be written 2lci6(x — x'). 
The expression for D is more complicated: one finds 

r, f > >\ ^i^^^'iV^y') fox X X \ 

Unm(x,x \ y,y ) = (2d„„i-d„+i,„i-()„_i,,„) . 

(4.8) 

The function R{x, x'; y, y') gives the overlap between the 
two directed intervals on the real line x — > cc' and y — s- 



y': for example, i?(l,5;4,9) = -i?(5,l;4,9) = 1. On 
substituting these expressions for C and D into Ea. H4.3|l . 
we obtain 



LHlN 



dxjdx'jdyjdy'{e 



1 1 ^ih{x — x') _j_ ^ — ih{x'-x') 



X j dtitG{x- x',t) J dt' it' G{y - y\t')e~^''-^'^/^-' 

X g^la-a'l/'cl J ^2R(x,x';y,y')/l^i _|_ ^-~2R{x,x' ; y ,y') /Ui _ ^ 



(4.9) 



Examining where the weight of the integrand lies with 
respect to the spatial integrals in Eq. I|4.9(l , one sees that 
the term in braces vanishes except in places where i? 7^ 0. 
We consider different types of contributions from these 
regions, and keep only those which are leading order for 
Lt ^ ^ei- First, consider regions in which jcc— y| ~ l^x but 
|x— x'l ^ lex- The small factor e.~\^~^ l/'oi is compensated 
by the first term in the braces \i\x' — y'\^ lc\. Then 



^-\x~x'\/l^i^-\y-y'\/l^i^2R{x,x'; y,y')/lai 



(4.10) 



g(-k-a|-|a:'-a'l)/'cl_ 

Since G{x,t) has a range in x of order Lt, the result- 
ing contribution to F{5B) is 0{Lt/L). Another con- 
tribution of the same order arises from regions where 
|x — y'l loi and \x' — y\ l^i- Subleading contribu- 
tions come from regions where all four spatial variables 
are within an elastic length of one another. These con- 
tributions are 0{lci/L). 

Keeping only the leading order terms, the expression 
for the correlation function has the much simplified form 



F{dB) 



NL 



rfa;(e*"^ + e"*''^) dttt dt' it 



X {G{x, t)G(x, t') + G{x, t)G{-x, t')) . (4.11) 
Using the symmetry of G{x, t) (see Eq. (|2.29|) ') one finds 



F{5B) 



9t 
NL 



dxe^'^^[f{x)Y 



where 



J{x) = -AttvI I dttlaiG{x,t) . 



(4.12) 



(4.13) 



B. Computing the correlation function 

In order to compare our theory for conductance fiuctu- 
ations with experiment, we need to be able to calculate 
F{6B) for various values of the temperature and param- 
eters vf and w. Although it is possible to use a computer 
to evaluate the form of F{B±) given in Eq. H4.12|l without 
further approximation, it is far easier to make progress 
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by calculating G{x, t) for a linearised dispersion relation. 
This approach is exact in the low-temperature regime de- 
fined by the condition cr{T) w cr(0), and we proceed to 
use it in our calculations. 

In the low temperature regime where the linearised dis- 
persion relation may be used, F(Bj_) has a scaling form. 
To make this apparent, it is helpful to recast equations in 
terms of dimensionless variables, characterising SB by b 
in place of h, and introducing x = x/ Lt and t ~ vpt/ L^. 
Writing G{x,t) = (27rLT)"^G(x, £) and /(Lxi) = f{x), 
for a linear dispersion relation, tjo{q,k) = qvFa{k), we 
have 



G{x, t) — exp 



-2a 



7r/a 

dk{l 



log \x — a{k)i\ — log 



— cos ak) 
7r[a(fc)i — x]/a{k) 



nTT/a 

X exp — ia dk{l 



sinh(T:[a{k)t — x]/a{k)) 
cos ak) sgn {x — a{k)i) > (4.14) 



and 



diilm{G{x,i)} . 



(4.15) 



Then the conductance autocorrelation function has the 
form 

9oLi 



F{5B) 
with scaling function 
C{b) = 



NL 



C {SB I Bo) 



dxe^'-[f{x)Y 



(4.16) 



(4.17) 



In this form F{SB) depends on temperature T and mag- 
netic field difference SB only through the scaling vari- 
ables Lt/L and 5B/Bq. The thermal length Lt plays 
the role of an inelastic scattering length, in the sense that 
it determines both the amplitude of conductance fluctu- 
ations and (through Bq) their correlation field. Such be- 
haviour is initally surprising, since Lj- is independent of 
interaction strength. In fact, of course, the form of the 
scaling function C (b) depends parametrically on interac- 
tion strength. 

For weak interactions this dependence of C{b) on k 
can be extracted analytically, as follows. First, note from 
Eg. j^in that a {k) = 1 + K/\k\. Also, in Eqs. f^H^ . 
(|4.15|) and 14.17|l . change variables from x,i to y,p with 
X — u/k and t — yp + y/n. Then 

\im G{y/K,p + y/K) = g{y,p) 

K.—>Q 



Try[p+l/k] 



exp < / dk{l — cos ak) 

I Jo 

X exp <ia dk{l — cos ak) sgn {y[p+l/k]) 




and 



0-5 8 B/Bg 1 



FIG. 6: C{SB / Bo) for narrow edge states and na — 0.6, 0.8, 
and 1. 



lim fiy/n) = f{y) 



TT 



dpp Im{g(y,_p)} 



The K-dependence of the scaling function is hence iso- 
lated for small k as 

C{b) = - / dyeM^yb/^)[fiy)? , (4.18) 

J -oo 

demonstrating that the amplitude of conductance fluc- 
tuations grows and that the correlation field shrinks as 
interactions are made weaker. In both cases, the varia- 
tion implies an inelastic scattering length that diverges 
as for weak interactions. Such a dependence of the 
inelastic scattering length on interaction strength is long- 
established in non-chiral, one-dimensional conductors. — 
In order to find the form of the scaling function and to 
study its K-dependence at general k, a three-dimensional 
numerical integration is necessary. We compute G(x,i), 
then /(i), and then the scahng function G(b) itself. 



C. Results 

We illustrate the form of the scaling function 
C{SB/Bq) for a range of parameter values in a sequence 
of three flgures. Its dependence on interaction strength 
Ka is shown for narrow edge states in Fig.|Hland for w = a 
in Fig. [71 In both cases, smaller interaction strength leads 
to a larger amplitude for conductance fluctuations and a 
smaller correlation field, as may be anticipated on the 
grounds that weaker interactions lead to a longer inelas- 
tic scattering length. In Fig.^G{SB / Bo) is shown for 
K = 50 and w — 4a, the parameter values suggested by 
the comparison of our conductivity calculations with ex- 
periment. We discuss experimental data on conductance 
fluctuations in Sec. IIVDI Finally, the increase in the 
amplitude of conductance fluctuations with dereasing k 
is illustrated in Fig. O 
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FIG. 8: C{5B/Bo) at w = 4a and Ka = 50. 




1 2 
log{K a) 

FIG. 9: Conductance fluctuation amplitude as a function of 
interaction strength at tt) = (full line), and asymptotic 
behaviour calculated analytically for small na (dashed line) 



D. Comparison with experiment and previous 
theory 

The exact treatment of disorder and interactions pro- 
vided by the calculations we have decribed presents an 
opportunity to test the standard theoretical treatment of 
conductance fluctuations, in which a single inelastic scat- 
tering length /in, or equivalently a scattering rate wp/^in 
is used as a cut-off in perturbation theory. For the chiral 
metal, such calculations have been described in Ref. 0. 
They yield a Lorcntzian scaling function 



F{dB) = 



NL 1 + 



(4.19) 



with z = 2Tr6Bliaa/^o. A comparison between the func- 
tional form we obtain for F{SB) and a Lorenztian is 
given in Fig. |S1 while the two functions are similar, the 
discrepancies are worth attention because they indicate 
behaviour which cannot be characterised by a single re- 
laxation time. A similar comparison can be made in 
the Fourier transformed domain, in terms of the func- 
tion f{x). To reproduce Eq. (|4.19() from our Eq. 14.12|l . 
we would require Un = Lt and 



(4.20) 



where exponential decay is indicative of a single lifetime 
hn/vp for excitations. The form we obtain for f{x) is 
shown in Fig. ^| The absence of a cusp at a; = indi- 
cates that there is of a range of relaxation times in the 
system. In addition, the fact that /(O) 7^ 1 is an interac- 
tion effect (from Eq. H3.10|l one sees that /(O) = a{0)/(Jo) 
not allowed for in the standard perturbative treatment. 

We close this section with a comparison between the 
experiments of Ref. and our results, using the same 
parameters, Ka = 50 and w = 4a, that provided a match 
for the behaviour of cr{T). For the experimental base 
temperature of T = 70mK, we use our approach to de- 
termine the amplitude of conductance fluctuations. As 
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a way to present the result, we then follow the experi- 
mental analysis!^ in using Eq. (|4.19|) to obtain a value 
for of 0.3/zm. The experimental value, extracted in 
the same way, is /jn ~ 1/im. Since the calculated ampli- 
tude of conductance fluctuations varies by several orders 
of magnitude over the range of parameter values we have 
investigated, and since no new adjustment of parameters 
was involved in our discussion of conductance fluctua- 
tions, we find the rough agreement between these two 
values of very encouraging. 

V. CONCLUSIONS 

In summary, for the system of weakly coupled quan- 
tum Hall edge states that we have studied, bosonisation 
provides a very complete treatment of the interplay be- 
tween electron-electron interactions and disorder. We 
have shown that interaction effects can account for the 
observed temperature dependence of interlayer conduc- 
tivity, provided we allow for finite edge state width and 



adopt a value for the edge state velocity that is rather 
smaller than previously supposed. We have investigated 
conductance fluctuations within the same theoretical ap- 
proach, showing how they are suppressed with increasing 
temperature, with a characteristic lengthscale cx T~^. 
Encouragingly, the same parameter values used to match 
the measured behaviour of conductivity reproduce ap- 
proximately the observed fluctuation amplitude. From a 
theoretical viewpoint, it is interesting that such dephas- 
ing effects can be generated from a description based on 
harmonic collective modes, simply via the nonlinear re- 
lation between boson and fermion operators. 
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